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Abstract.  Phlebotomus  papatasi  sand  flies  are  among  the  primary  vectors  of  Leishmania  major  parasites  from 
Morocco  to  the  Indian  subcontinent  and  from  southern  Europe  to  central  and  eastern  Africa.  Antibody-based  immunity 
to  sand  fly  salivary  gland  proteins  in  human  populations  remains  a  complex  contextual  problem  that  is  not  yet  fully 
understood.  We  profiled  the  immunoreactivities  of  plasma  antibodies  to  sand  fly  salivary  gland  sonicates  (SGSs)  from 
229  human  blood  donors  residing  in  different  regions  of  sand  fly  endemicity  throughout  Jordan  and  Egypt  as  well  as 
69  US  military  personnel,  who  were  differentially  exposed  to  P.  papatasi  bites  and  L.  major  infections  in  Iraq.  Compared 
with  plasma  from  control  region  donors,  antibodies  were  significantly  immunoreactive  to  five  salivary  proteins  (12,  26, 

30,  38,  and  44  kDa)  among  Jordanian  and  Egyptian  donors,  with  immunoglobulin  G4  being  the  dominant  anti-SGS 
isotype.  US  personnel  were  significantly  immunoreactive  to  only  two  salivary  proteins  (38  and  14  kDa).  Using  k-means 
clustering,  donors  were  segregated  into  four  clusters  distinguished  by  unique  immunoreactivity  profiles  to  varying 
combinations  of  the  significantly  immunogenic  salivary  proteins.  SGS-induced  cellular  proliferation  was  diminished 
among  donors  residing  in  sand  fly-endemic  regions.  These  data  provide  a  clearer  picture  of  human  immune  responses  to 
sand  fly  vector  salivary  constituents. 


INTRODUCTION 

The  causative  agents  of  leishmaniasis  are  protozoan  para¬ 
sites  of  the  genus  Leishmania.  All  species  within  the  genus 
share  a  common  lifecycle  that  involves  an  extracellular  infec¬ 
tive  stage  present  in  the  sand  fly  vector  followed  by  an  intra¬ 
cellular  stage  of  replication  within  a  vertebrate  host.  Estimates 
of  annual  global  disease  incidence  are  approximately  0.2- 
0.4  and  0.7-1 .2  million  cases  of  active  visceral  and  cutaneous 
leishmaniasis,  respectively.1  Leishmaniasis  is  primarily  found  in 
tropical  and  subtropical  regions,  where  vector-competent  sand 
fly  species  of  the  family  Psychodidae  reside.  Medically  relevant 
genera  of  the  sand  fly  subfamily  Phlebotominae  are  distin¬ 
guished  by  distinct  geographical  regions  of  endemicity: 
Phlebotomus  in  the  Old  World  (the  Mediterranean,  Middle 
East,  and  Africa)  and  Lutzomyia  in  the  New  World  (Central 
and  South  America).2 

Currently,  no  efficacious  vaccine  exists  to  combat  infection 
or  transmission  of  Leishmania.  In  lieu  of  a  vaccine  directed  at 
the  parasites,  attention  has  focused  on  another  aspect  of  the 
disease  cycle:  the  sand  fly  vector.  One  hypothesis  is  that, 
through  pre-exposure  to  sand  fly  vector  salivary  components, 
individuals  are  immunized  to  establish  a  robust  immunologi¬ 
cal  memory  against  vector  salivary  factors.  This  response  may 
generate  a  hostile  microenvironment  for  parasite  transmission 
on  reexposure.  For  every  species  of  Leishmania  parasite,  only 
a  few  sand  fly  species  are  capable  of  vectoring  a  specific 
Leishmania  species,3  resulting  in  a  very  restrictive  and  unique 
vector-parasite  exposure  profile.  Key  to  the  downstream 
development  of  a  vaccine  is  the  careful  examination  of  host 
immunological  responses  to  vector  salivary  proteins  among 
persons  naturally  pre-exposed  to  sand  fly  bites.4-5 
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Many  studies  have  sought  to  investigate  the  serological 
markers  of  hosts,  both  animal  and  human,  against  sand  fly 
species,6-14  although  the  problem  remains  a  highly  complex 
and  contextual  matter  warranting  deeper  novel  analysis.  To 
better  characterize  acquired  human  host  immune  responses  to 
sand  fly  saliva,  we  collected  peripheral  blood  samples  from 
229  people  who  reside  in  Middle  Eastern  regions  (Egypt  and 
Jordan)  of  variable  endemic  disease  ecologies  for  cutaneous 
leishmaniasis  and  its  predominant  Old  World  vector,  the  sand 
fly  P.  papatasi.  Additionally,  samples  were  collected  from 
69  US  military  personnel,  differentially  exposed  to  P.  papatasi 
bites  and  L.  major  infections,  to  aid  in  comparisons  with  immu- 
nologically  naive  individuals.  Our  work  sought  to  understand 
host  population  variation  with  reference  to  antivector  saliva 
immunity  and  identify  salivary  proteins  for  later  investigation 
as  potential  exposure  markers. 

METHODS 

Study  sites  and  human  peripheral  blood  donors.  Six  sites  in 
the  Middle  East  and  North  Africa  (MENA)  were  chosen  for 
collection  of  229  human  peripheral  blood  samples  during  late 
August  and  early  September  of  2007.  Table  1  shows  regional 
donor  profiles.  Geospatial  coordinates  for  these  regions  were 
as  previously  described.15  Each  site  was  selected  based  on  eco- 
logical  characteristics  in  Jordan  -  and  Egypt-  -  pertaining 
to  active  cases  of  leishmaniasis  transmission  and  the  presence 
of  P.  papatasi  vectors.  These  sites  included  Cairo  (CA),  Egypt 
and  Amman  (AM),  Jordan  (a  lack  of  sand  fly  vectors);  Aswan 
(AW),  Egypt  and  Malka  (MA),  Jordan  (regions  of  endemic  P. 
papatasi  populations  but  that  lack  clinical  reports  of  L.  major 
infections);  and  North  Sinai  (NS),  Egypt  and  Swaymeh  (SW), 
Jordan  (regions  of  endemic  P.  papatasi  population  as  well  as 
regions  with  reports  of  active  L.  major  transmission).  Blood 
samples  from  regional  Egyptian  and  Jordanian  donors  were 
drawn  from  healthy  adults,  none  of  whom  presented  with 
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Table  1 


Summary  of  human  peripheral  blood  donors 


Abbreviation 

Country 

P.  papatasi  endemicity* 

L.  major  transmission  reports* 

Number  of  blood  donors 

MENA  sitesf 
Cairo 

CA 

Egypt 

37 

Aswan 

AW 

Egypt 

+ 

- 

40 

North  Sinai 

NS 

Egypt 

+ 

+ 

34 

Amman 

AM 

Jordan 

- 

- 

37 

Malka 

MA 

Jordan 

+ 

- 

41 

Swaymeh 

SW 

Jordan 

+ 

+ 

40 

US  military!: 
Control 

United  States 

P.  papatasi  bite  exposure  § 

L.  major  infection§ 

29 

Bitten 

United  States 

+ 

- 

17 

Infected 

United  States 

+ 

+ 

23 

*  Marks  indicate  the  presence  (+)  or  absence  (-)  of  endemic  P.  papatasi  populations  and  known  cases  of  active  L.  major  transmission  in  the  region. 
tMENA  locations  in  which  human  peripheral  blood  donors  resided. 

tThe  condition  of  US  military  personnel  stationed  in  the  United  States  (control)  or  deployed  in  Iraq  (bitten  and  infected)  during  2007. 

§  Marks  indicate  yes  (+)  or  no  (-)  as  to  whether  US  military  personnel  reported  exposure  to  sand  fly  bites  or  there  was  a  confirmed  diagnosis  of  cutaneous  leishmaniasis  while  deployed 
in  Iraq. 


active  lesions.  Peripheral  blood  samples  were  similarly  col¬ 
lected  from  69  US  military  personnel  at  Walter  Reed  Army 
Medical  Center,  Washington,  DC.  US  donors  were  interviewed 
for  histories  of  suspected  sand  fly  bite  exposure  or  presented 
with  clinically  confirmed  cases  of  cutaneous  leishmaniasis  but  were 
untreated  before  blood  sampling.  All  donors  were  >  18  years 
of  age,  and  the  number  of  male  and  female  participants  was 
approximately  equal.  All  donors  signed  documents  of  informed 
consent  per  local  and  US  laws,  and  samples  were  stripped  of  all 
personal  identifying  information  before  processing. 

Human  host  blood  collection.  Ten  to  fifteen  milliliters 
blood  were  collected  in  yellow  top  Vacutainer  blood  collec¬ 
tion  tubes  with  ACD  solution  A  anticoagulant  (BD  Biosci¬ 
ences,  San  Jose,  CA).  Tubes  were  kept  at  4°C  and  centrifuged 
at  3,000  rpm  for  10  minutes.  The  aqueous  top  plasma  layer 
was  removed  and  stored  at  -20°C.  Remaining  mixed  leuko¬ 
cytes  and  red  blood  cells  were  segregated  using  Lymphocyte 
Separation  Media  (Cellgro,  Manassas,  VA)  following  the  man¬ 
ufacturer’s  protocol.  Remaining  cells  were  suspended  in  red 
blood  cell  lysis  buffer  (0.15  M  ammonium  chloride,  10  mM 
potassium  bicarbonate,  and  0.1  mM  (ethylenedinitrilo) 
tetraacetic  acid  [EDTA]),  incubated  for  5  minutes  on  ice, 
washed  with  lx  phosphate  buffered  saline  (PBS;  137  mM 
NaCl,  2.7  mM  KC1,  10  mM  Na2HP04,  and  1.8  mM  KH2P04), 
and  resuspended  in  standard  RPMI-1640  media  (Cellgro)  sup¬ 
plemented  with  10%  fetal  bovine  serum  (FBS),  100  U/mL 
penicillin,  100  pg/mL  streptomycin,  and  2  mM  L-glutamine 
(referred  to  as  RPMI-C).  Mixed  cells  in  RPMI-C  media  were 
combined  1:1  in  freezing  media  (20%  dimethyl  sulfoxide 
[DMSO],  20%  FBS,  and  60%  RPMI-C),  slow  frozen  overnight 
at  -80  °C  in  an  isopropanol  chamber,  transferred  to  dry  ice  for 
transport,  and  ultimately  transferred  to  liquid  nitrogen  for 
long-term  storage  before  use  in  cell  proliferative  assays. 

Sand  fly  colony.  P.  papatasi  sand  flies  used  in  this  study 
were  derived  from  the  Israeli  strain  maintained  at  the  Labo¬ 
ratory  of  Malaria  and  Vector  Research,  National  Institutes  of 
Health  (Rockville,  MD)  and  the  University  of  Notre  Dame 
(Notre  Dame,  IN)  as  previously  described.15 

Salivary  gland  sonicates.  Salivary  glands  from  1-  to  2-week- 
old  unfed  female  P.  papatasi  sand  flies  were  dissected  intact 
under  1  x  PBS  and  flash  frozen  at  -80°C  for  storage.  To  pro¬ 
cure  salivary  gland  protein  sonicates  (SGSs),  salivary  glands 
were  subjected  to  three  cycles  of  rapid  freeze-thaw  (-80°C  for 
2  minutes  and  then  a  37°C  water  bath  for  1  minute),  centrifuged 
for  1  minute  at  14,000  rpm,  sonicated  in  a  room  temperature 


water  bath  for  1  minute  at  80%  power  using  an  Ultrasonic 
Processor  (Cole  Parmer,  Vernon  Hills,  IL),  and  then,  filter- 
sterilized  by  centrifugation  at  14,000  rpm  for  1  minute  through 
a  0.22-pm  Ultrafree  centrifugal  filter  unit  (Millipore,  Billerica, 
MA).  Total  protein  quantities  were  assessed  by  standard  curve 
to  bovine  serum  albumin  dilutions  using  the  Protein  A280 
feature  of  a  NanoDrop  1000  microvolume  UV-Vis  spectropho¬ 
tometer  (Thermo  Scientific,  Waltham,  MA).  A  whole-salivary 
gland  equivalent  concentration  was  considered  approximately 
1.5  pg  based  on  estimates  of  complete  lobe  pairs  from  sugar-fed 
female  Lu.  longipalpis  salivary  gland  protein  concentrations.30 

Mass  spectrometry  and  protein  sequence  alignments. 
Twenty  salivary  gland  protein  equivalents  from  P.  papatasi 
SGS  in  lx  PBS  were  boiled  at  85°C  for  5  minutes  with 
Tris-glycine  sodium  dodecyl  sulfate  (SDS)  sample  buffer  and 
NuPAGE  reducing  agent  (Novex,  Carlsbad,  CA).  Samples 
were  applied  to  a  pre-cast  4-20%  stacking  Tris-glycine  gel 
(Novex)  with  a  PageRuler  pre-stained  protein  ladder 
(Fermentas,  Pittsburgh,  PA)  on  the  side  and  electrophoresed 
at  125  V  for  90  minutes.  Salivary  proteins  were  run  in  single 
dimension  electrophoresis  to  allow  cross-comparison  with 
immunoblots.  The  gel  was  bathed  in  fixation  solution  (50% 
methanol  and  10%  acetic  acid)  overnight  and  then  stained  for 
4  hours  in  1%  Coomassie  Blue  R-250  (Bio- Rad,  Hercules, 
CA).  The  gel  was  destained  overnight  in  destaining  solution 
(5%  methanol  and  7.5%  acetic  acid)  and  stored  in  7%  acetic 
acid  for  imaging  on  a  gel  doc  Universal  Hood  II  (Bio-Rad, 
Hercules,  CA).  Ten  intensely  stained  protein  gel  bands  were 
cut  precisely  using  clean,  sterile  scalpels.  Gel  slices  were  pre¬ 
pared  by  trypsin  digestion  for  mass  spectrometry  (MS)  as  pre¬ 
viously  described.31  Samples  were  desalted  by  Ziptip 
(Millipore,  Billerica,  MA)  before  liquid  chromatography  (LC) 
/MS/MS  analysis  performed  on  a  5500QTrap  operating  in  ion- 
trap  mode  (ABSciex,  Farmingham,  MA)  connected  to  an 
Eksigent  (AB  Sciex)  Nano  Ultra  UHPLC  system.  Samples  were 
acquired  as  described  32  Briefly,  each  gel  band  was  analyzed  in 
technical  duplicate  and  separated  over  a  gradient  from  2%  to 
35%  acetonitrile  (0.1%  formic  acid)  on  a  100-pm  x  100-mm 
C18BEH  column  (Waters,  Milford,  MA)  running  at  600  nL/ 
min  for  1  hour.  A  TOP8  data-dependent  method  was  used  for 
spectral  acquisition.  Peak  lists  were  subjected  to  database 
search  against  all  available  Phlebotomus  and  Lutzomyia  spp. 
protein  sequences  available  through  the  US  National  Center 
for  Biotechnology  Information  (NCBI)33  combined  with  a  list 
of  common  contaminants  (e.g.,  trypsin,  keratin,  etc.).  False 
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discovery  rates  (FDRs)  were  determined  by  the  target-decoy 
approach.34  The  Paragon  algorithm  within  the  Protein  Pilot 
software  suite  (AB  Sciex,  Farmingham,  Massachusetts,  USA) 
was  used  for  database  searches,32  and  proteins  within  the 
95%  FDR  interval  were  chosen  for  additional  analysis.  FDR 
descriptions  for  combined  protein,  peptide,  and  spectral  level 
analyses  are  included  in  Supplemental  Figure  1.  Alignment 
analyses  of  published  protein  sequences  were  conducted 
using  Clustal  Omega35  or  LAlign36  web-based  programs  with 
default  settings. 

Inununoblots.  Western  immunoblots  to  examine  human 
plasma  donor  antibodies  reactive  against  SGS  were  conducted 
using  a  modified  version  of  a  previously  reported  protocol.37 
Briefly,  30  SGSs  were  boiled  and  separated  on  pre-cast 
10%  Tris-glycine  SDS-polyacrylamide  gel  electrophoresis 
(PAGE)  gels  (Novex)  with  a  PageRuler  pre-stained  protein 
ladder  (Fermentas)  on  the  side  and  transferred  to  nitrocellu¬ 
lose  membranes  using  an  XCell  II  blot  module  (Novex). 
Membranes  were  blocked  for  1  hour  at  room  temperature 
(RT)  in  blocking  buffer  (0.1%  Tween  20/lx  tris-buffered 
saline  and  5%  dried  milk),  washed  three  times  in  1  x  TBS, 
and  set  in  a  Mini-Protean  II  multiscreen  apparatus  (Bio- 
Rad).  Blood  plasma  from  human  donors  was  diluted  1:100  in 
blocking  buffer  and  exposed  to  the  membranes.  A  single  com¬ 
mon  donor  was  exposed  to  Israeli  colony  P.  papatasi  bites  24 
hours  before  blood  collection,  which  served  as  a  normalizing 
control  on  each  blot.  Membranes  and  plasma  were  incubated 
overnight  at  4°C,  washed  three  times  in  lx  TBS,  and  then, 
exposed  to  a  secondary  alkaline  phosphatase-conjugated  goat 
anti-human  Fc  antibody  (Zymed,  South  San  Francisco,  CA) 
diluted  1:5,000  in  blocking  buffer.  Membranes  were  devel¬ 
oped  with  colorimetric  substrate  One  Step  NPT-BCIP 
(Thermo  Scientific).  Membranes  were  imaged  using  a  gel  doc 
Universal  Hood  II  (Bio-Rad).  White  balance  for  all  images 
was  normalized  to  a  single  immunoblot  (Figure  1A)  using 


Photoshop  CS5  software  (Adobe,  San  Jose,  CA)  to  equalize 
background  intensity.  Because  of  non-specific  colorimetric 
artifacts  in  the  lanes  of  many  immunoblots,  the  use  of  auto¬ 
mated  image  analysis  software,  such  as  ImageJ,  to  score  the 
plasma  antibody  reactivity  to  SGS  protein  bands  proved 
impossible,  because  a  universal  baseline  was  not  applicable 
to  all  lanes,  even  in  a  single  immunoblot.  Three  people  man¬ 
ually  scored  each  blot,  and  they  were  independent  observers 
blinded  to  the  regional  identities  of  the  plasma  donors  applied 
to  each  immunoblot.  Observers  cataloged  the  banding  pattern 
and  intensity  of  each  blot  based  on  a  numeric  scale  (0  =  no 
bands  at  a  given  kilodalton  value;  1  =  low-intensity  band;  2  = 
moderate-intensity  band;  3  =  high-intensity  band)  relative  to 
the  positive  control  plasma  lane  appearance.  The  scores  were 
combined  based  on  a  uniform  rule  set  favoring  the  majority 
opinion  at  each  kilodalton  range.  The  combined  score  calcu¬ 
lations  were  checked  and  verified  by  each  investigator  inde¬ 
pendently,  providing  the  most  unbiased  dataset  possible. 

Enzyme-Linked  Immunosorbent  Assay.  Isolated  SGSs  in 
1  x  PBS  were  applied  at  1.5  pg  per  well  onto  a  MaxiSorp  flat- 
bottom  96-well  plate  (Nunc,  Penfield,  NY)  and  incubated 
overnight  at  4°C.  Some  wells  on  each  plate  were  coated 
instead  with  purified  human  immunoglobulin  G  (IgG)  or  IgE 
(Abeam,  Cambridge,  MA)  at  known  concentrations  for  1:100, 
1:1,000,  and  1:10,000  dilutions,  which  served  as  standards  for 
calculations  of  human  donor  anti-SGS  antibody  concentra¬ 
tions.  IgG  standards  (Abeam)  were  composed  of  a  mixture  of 
all  subclasses  (IgGl,  -2,  -3,  and  -4)  at  known  percentages  of 
the  total.  Plates  were  washed  with  0.05%  Tween-20/1  x  PBS 
and  blocked  with  5%  bovine  serum  albumin  (BSA)  in  1  x  PBS 
for  1  hour  at  RT.  Human  plasma  in  blocking  solution  was 
applied  at  a  1:100  dilution  and  incubated  overnight  at  4°C. 
Mouse  secondary  biotinylated  anti-human  antibodies,  specific 
for  IgGl,  -2,  -3,  and  -4  and  IgE  isotypes  (Sigma-Aldrich, 
Saint  Louis,  MO),  were  applied  at  a  1:100  dilution  in  blocking 
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Figure  1.  Characterization  of  whole  P.  papatasi  salivary  gland  protein  extract  sonicates  and  human  host  antibody  reactions.  (A)  Immunoblot 
examples  of  plasma  donors  from  three  regions  in  Egypt  and  differentially  exposed  US  military  personnel  with  antibody  specificity  to  P.  papatasi 
salivary  proteins.  Lanes  are  designated  by  regional  plasma  donor  codes  (CA,  AW,  and  NS)  or  differential  exposure  codes  (Con  =  Control,  Bit  = 
Bitten,  and  Inf  =  Infected).  (B)  Comparison  of  immunoblot  example  and  whole  P.  papatasi  salivary  gland  protein  sonicate.  Lane  2  displays  an 
immunoblot  image  of  a  North  Sinai  plasma  donor  with  broad  antibody  specificities  to  P.  papatasi  salivary  gland  proteins.  Lane  3  displays  a 
Coomassie  Brilliant  Blue  R-250-stained  P.  papatasi  whole-salivary  gland  protein  sonicate  run  on  a  denaturing  SDS-PAGE  gel.  Lanes  1  and  4  are 
PageRuler  protein  MW  ladders  (Fermentas). 
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buffer,  covered  from  light,  and  incubated  overnight  at  4°C. 
Plates  were  incubated  with  streptavidin-conjugated  alkaline 
phosphatase  (Sigma- Aldrich)  diluted  1:500  in  TBS-MgCl2/l% 
BSA  solution  for  2  hours  at  RT  and  covered  from  light.  After  a 
final  wash,  para-nitrophenyl  phosphate  (pNPP)  substrate  buffer 
(Sigma- Aldrich)  was  applied,  and  the  plates  were  incubated  for 
30  minutes  at  RT  and  read  at  405  nm  on  a  SpectraMax  M2 
microplate  reader  (Molecular  Devices,  Sunnyvale,  CA).  IgG 
and  IgE  standards  were  used  to  calculate  donor  anti-SGS  anti¬ 
body  concentrations  by  standard  curve  fitting. 

Proliferation  assays.  Monocytes  from  cryopreserved  mixed 
leukocyte  samples  of  human  donors  were  collected  by  adher¬ 
ence  to  six-well  cell  culture  plates  (Costar,  Manassas,  VA) 
after  quick  thaws  (1  minute  in  a  37  °C  water  bath),  washing 
in  RPMI-C  media,  and  a  2-hour  incubation  at  37°C  under 
5%  C02.  Cells  were  stimulated  every  other  day  for  1  week 
with  recombinant  human  granulocyte/macrophage  colony- 
stimulating  factor  (GM-CSF;  400  U  per  1  x  106  monocytes) 
and  IL-4  (160  U  per  1  x  106  monocytes)  to  generate  monocyte- 
derived  dendritic  cells  (MDDCs).  Cells  were  harvested  by 
repeated  washing  with  0.5  mM  EDTA  in  lx  PBS,  replated  in 
a  96-well  flat-bottom  cell  culture  plate  at  a  concentration  of 
4  x  104  MDDCs  per  well  in  RPMI-C  media,  and  then,  allowed 
to  rest  overnight  at  37  °C  under  5%  C02.  Cells  were  stimulated 
with  SGS  overnight  at  a  concentration  of  one  salivary  gland 
equivalent  (1.5  pg)  per  well.  Negative  control  wells  received 
only  RPMI-C  media,  whereas  positive  control  wells  received 
400  ng  phytohemagglutinin  (PHA;  Sigma-Aldrich)  per  well. 
Mixed  leukocytes  from  the  same  donors  were  coincubated 
with  stimulated  MDDCs  at  a  concentration  of  8  x  104  per  well 
at  37°C  under  5%  C02  for  4  days;  1  pCi  tritiated  thymidine 
(Perkin-Elmer,  Waltham,  MA)  in  RPMI-C  was  added  to  each 
well  for  the  final  24  hours  of  coculture.  Mixed  cultures  were 
harvested  using  a  scintillation  harvester  and  filter  mat  paper 
(Brandel,  Gaithersburg,  MD)  per  the  manufacturer’s  proto¬ 
col.  Scintillation  fluid  allowed  p-particle  decay  readings  on  a 
Wallac  MicroBeta  TriLux  microplate  liquid  scintillation 
counter  (Perkin-Elmer).  Final  values  are  reported  as  the 
stimulation  index— the  ratio  of  p-counts  (3H-thymidine 
uptake)  between  unstimulated  and  SGS  pre-stimulated 
MDDCs  and  mixed  leukocyte  cocultures. 

Luminex  cytokine  quantification.  Culture  supernatants 
from  the  proliferation  reactions  were  collected  just  before 
3H-thymidine  addition.  Supernatants  were  prepared  for  quan¬ 
titative  cytokine  and  chemokine  fluorescence-based  analysis 
using  a  Human  Cytokine  LINCOplex  Premixed  96-well  kit 
(Millipore,  Billerica,  MA)  following  the  manufacturer’s  pro¬ 
tocols.  Plates  were  read  on  a  Luminex  200  machine,  and  the 
mean  fluorescence  intensity  (MFI)  was  calculated  following 
the  manufacturer’s  protocols  for  internal  normalization.  Final 
MFI  values  for  each  of  21  cytokines  or  chemokines  in  super¬ 
natants  from  each  culture  exposed  in  vitro  to  SGS  were  sub¬ 
sequently  normalized  to  supernatant  values  of  paired  negative 
control  cultures  left  unexposed  to  SGS. 

Clustering  analysis.  K-means  clustering  was  performed 
using  Cluster  v.3.0.38  Data  were  formatted  as  sample  (plasma 
donor)  and  gene  (SGS  protein  kilodaltons)  tables  to  generate 
k  =  4  clusters  using  dataset  means  measured  by  Euclidean 
distance.  Visualization  and  data  export  were  achieved  with 
Java  Treeview  v. 1.1.4. 39 

Statistics.  Using  GraphPad  Prism  v.5.0  (Graphpad  Software, 
San  Diego,  CA),  combined  immunoblot  scoring  data  were  ana¬ 


lyzed  by  one-way  analysis  of  variance  (ANOVA)  and 
Bonferroni  multiple  comparison  post-tests  performed  at  each 
SGS  protein  kilodalton  for  every  regional  plasma  donor 
group  at  a  95%  confidence  interval.  A  simple  calculation  was 
used  to  describe  immunogenicity  of  each  salivary  protein: 
(probability  of  donor  plasma  antibody  specificity  to  a  particu¬ 
lar  salivary  protein  within  a  regional  sample  group)  x  (aver¬ 
age  band  intensity  among  that  regional  group)  =  arbitrary 
units  (a.u.)  of  immunogenicity  between  zero  and  two.  The 
same  software  was  used  to  calculate  Pearson  correlation  by 
two-tailed  distribution  using  percentages  of  individuals  from 
each  single  region  that  clustered  into  the  four  k-means  clusters 
as  well  as  Student’s  t  test  comparing  regional  groups  mixed 
leukocyte  proliferation  data. 

RESULTS 

Host  antibody  reactivity  to  vector  salivary  proteins.  To  better 
understand  human  anti-sand  fly  salivary  antibody  presence 
based  on  differential  exposure  to  wild  P.  papatasi  sand  flies, 
six  sites  were  selected  as  collection  points  for  human  host 
peripheral  blood  samples.  Each  regional  site,  three  sites 
in  Egypt  and  three  sites  in  Jordan,  was  categorized  into 
different  P.  papatasi  and  leishmaniasis  eco topes  (Table  1). 
Additionally,  US  military  personnel  represented  responses 
by  individuals  immunologically  naive  to  newly  encountered 
sand  fly  saliva.  Plasma  from  a  total  of  111  Egyptian  or 
118  Jordanian  residents  and  69  US  military  personnel  was  tested 
against  PAGE-separated  SGS  from  inbred  female  P.  papatasi 
by  Western  immunoblot  (Figure  1A  and  Supplemental  Figure 
2).  Concurrently,  PAGE-separated  SGS  revealed  11  heavily 
expressed  salivary  proteins  (~12,  14,  26,  28,  30,  38,  40,  42,  44, 
46,  and  58  kDa),  with  1  modestly  expressed  protein  of 
~68  kDa  (Figure  IB).  Human  donor  blood  plasma  antibodies 
displayed  differential  specificities  for  those  proteins  (Figure  1). 
Images  of  immunoblots  were  scored  by  signal  intensity  of 
plasma  antibodies  specific  for  P.  papatasi  salivary  proteins  at 
each  approximate  molecular  weight  (MW).  All  final  immuno¬ 
blot  scores  at  each  MW  for  Egyptian  and  Jordanian  donors 
were  compared  with  those  scores  of  donors  residing  in  control 
regions,  where  no  sand  fly  vectors  were  reported  to  be  present. 
US  military  donor  scores  were  compared  with  those  scores  of 
control  donors  who  denied  any  travel  to  P.  papatasi-e ndemic 
regions.  Antibody  specificities  for  five  salivary  proteins  (~12, 
26,  30,  38,  and  44  kDa)  were  found  to  be  significantly  dif¬ 
ferent  from  control  regions  within  one  or  both  Middle  East¬ 
ern  nations,  whereas  antibody  specificity  was  significantly 
different  compared  with  the  control  group  for  only  two  salivary 
proteins  (~14  and  38  kDa)  among  Iraq-based  US  donors 
infected  with  Leishmania  (P  <  0.05)  (Table  2).  Notably,  among 
Middle  Eastern  donors,  there  seemed  to  be  very  little,  if  any, 
specificity  for  salivary  proteins  of  ~15  or  21  kDa  (data  not 
shown)  and  some  minor  differential  recognition  of  an  ~36  kDa 
protein  that  was  not  observed  to  be  a  majorly  expressed  pro¬ 
tein  in  the  P.  papatasi  laboratory  colony  SGS  (Figure  IB). 
However,  a  salivary  protein  at  ~14  kDa  was  significantly  differ¬ 
entially  recognized  by  antibodies  from  US  military  donors. 

Antibody  specificity  for  salivary  proteins  can  be  summa¬ 
rized  using  a  novel  calculation  of  immunogenicity  (Figure  2). 
Regarding  Egyptian  and  Jordanian  donors,  plasma  displayed 
the  greatest  immunogenicity  against  salivary  proteins  at  ~12, 
30,  and  38  kDa,  regardless  of  geographic  origin.  Greater  relative 
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Table  2 


Characteristics  of  salivary  proteins  with  significant  differences  in  host  antibody  specificity  profiles 


Salivary  proteins  (estimated  MW) 

12  kDa 

26  kDa 

30  kDa 

38  kDa 

44  kDa 

MENA* 

Percentage  positive  seraf 

Egypt 

CA  (-/-) 

43.24 

27.03 

37.84 

45.95 

8.11 

AW  (+/-) 

90.00f 

70.00* 

90.00* 

95.00* 

42.50* 

NS  (+/+) 

79.41* 

55.88* 

88.24* 

79.41* 

41.18* 

Jordan 

AM  (-/-) 

56.76 

21.62 

59.46 

70.27 

13.51 

MA  (+/-) 

58.54 

21.95 

78.05* 

73.17 

12.20 

SW  (+/+) 

77.50f 

35.00 

90.00* 

60.00 

45.00* 

Average  blot  intensity§ 

Egypt 

CA  (-/-) 

0.51  ±  0.08 

0.38  ±  0.13 

0.62  ±  0.15 

0.59  ±  0.12 

0.14  ±0.11 

AW  (+/-) 

1.31  ±0.15* 

0.90  ±0.11* 

1.83  ±  0.16* 

1.33  ±0.11* 

0.70  ±  0.12* 

NS  (+/+) 

1.26  ±0.18* 

0.96  ±  0.17* 

2.09  ±  0.18* 

1.38  ±  0.17* 

0.84  ±  0.15* 

Jordan 

AM  (-/-) 

0.66  ±  0.10 

0.24  ±0.14 

0.86  ±  0.15 

0.97  ±  0.08 

0.22  ±0.11 

MA  (+/-) 

0.76  ±  0.06 

0.29  ±  0.15 

1.37  ±  0.16* 

1.15  ±0.10 

0.13  ±0.12 

SW  (+/+) 

1.08  ±0.14* 

0.41  ±  0.15 

1.53  ±0.15* 

0.94  ±  0.10 

0.75  ±  0.12* 

Immunogenicity  (a.u.)1 

Egypt 

CA  (-/-) 

0.22 

0.10 

0.24 

0.27 

0.01 

AW  (+/-) 

1.18* 

0.63* 

1.64* 

1.26* 

0.30* 

NS  (+/+) 

1.00* 

0.53* 

1.84* 

1.10* 

0.35* 

Jordan 

AM  (-/-) 

0.38 

0.05 

0.51 

0.68 

0.03 

MA  (+/-) 

0.44 

0.06 

1.07* 

0.84 

0.02 

SW  (+/+) 

0.83* 

0.14 

1.37* 

0.56 

0.34* 

12  kDa 

14  kDa 

26  kDa 

30  kDa 

38  kDa 

US  military* 

Percentage  positive  seraf 

United  States 

Control 

48.28 

51.72 

44.83 

48.28 

48.28 

Bitten 

23.53 

52.94 

64.71 

29.41 

76.47 

Infected 

34.78 

13.04* 

30.43 

21.74 

73.91* 

Average  blot  intensity§ 

United  States 

Control 

0.52  ±  0.57 

0.53  ±  0.53 

0.50  ±  0.60 

0.64  ±  0.77 

0.48  ±  0.51 

Bitten 

0.24  ±  0.44 

0.71  ±  0.85 

0.74  ±  0.62 

0.29  ±  0.47 

1.00  ±  0.71 

Infected 

0.52  ±  0.85 

0.30  ±  0.82* 

0.48  ±  0.79 

0.30  ±  0.70 

1.70  ±  1.19* 

Immunogenicity  (arbitrary  units)1]! 

United  States 

Control 

0.25 

0.28 

0.22 

0.31 

0.23 

Bitten 

0.06 

0.37 

0.48 

0.09 

0.76 

Infected 

0.18 

0.04* 

0.15 

0.07 

1.25* 

*  Collection  site  abbreviations  and  conditional  marks  correspond  to  those  abbreviations  and  marks  presented  in  Table  1  (e.g.,  +/-  refers  to  a  region  with  known  endemic  P.  papatasi  populations 
but  no  reports  of  L.  major  infections  or  for  US  military  personnel,  donor  reported  history  of  sand  fly  bites  but  there  is  no  confirmed  diagnoses  of  cutaneous  leishmaniasis), 
t  Percentage  of  serum  donors  from  the  given  location  with  antibody  recognition  of  salivary  proteins  at  the  estimated  MW. 

t Those  regions  with  significant  (P  <  0.05)  differences  in  antibody  reactivity  compared  with  national  negative  control  regions  or  US  personnel  who  were  not  bitten  by  sand  flies. 

§  Western  immunoblots  of  human  sera  against  P.  papatasi  SGS  were  scored  manually  by  three  independent  observers  using  a  scale:  0  =  no  reaction;  1  =  weak  reaction;  2  =  moderate  reaction; 
3  =  intense  reaction  (all  judged  relative  to  a  constant  positive  control). 

1  Immunogenicity  here  is  calculated  as  (probability  of  donor  serum  antibody  specificity  to  a  particular  salivary  protein  within  a  regional  sample  group)  x  (average  band  intensity  among  that 
regional  group)  =  a.u.  of  immunogenicity. 


difference  to  control  regions  was  more  readily  apparent  among 
plasma  from  Egyptian  donors.  Jordanian  donor  plasma  from 
Amman,  a  control  region,  exhibited  much  higher  immunoge¬ 
nicity  to  those  three  salivary  proteins  compared  with  donor 
plasma  from  Cairo,  the  Egyptian  control  region  (Figure  2  and 
Table  1).  This  increased  antibody  specificity  for  three  highly 
immunogenic  salivary  proteins,  in  part,  accounts  for  the  lack  of 
statistical  significance  in  testing  differences  between  intrare- 
gional  profiles  for  Jordanian  donors.  The  percentage  of  indi¬ 
viduals  reacting  to  the  five  significant  salivary  proteins  was 
greater  among  individuals  residing  in  AW  compared  with  MA 
(Table  2).  However,  between  NS  and  SW,  the  percent  of  per¬ 
sons  reacting  to  salivary  proteins  of  ~12,  30,  and  44  kDa  was 


similar,  although  more  Egyptian  donors  displayed  antibodies 
against  salivary  proteins  of  ~26  and  38  kDa. 

Considering  the  immunogenicity  of  salivary  proteins  to  US 
military  donors,  values  for  the  ~38  kDa  salivary  protein  were 
substantially  greater  than  those  values  of  all  other  proteins 
(Figure  2).  Although  64%  of  bite-exposed  individuals  displayed 
specificity  for  the  -26  kDa  salivary  protein  (Table  2),  signal 
intensity  and  immunogenicity  were  not  significantly  different 
from  signal  intensity  and  immunogenicity  of  control  donor 
plasma.  Although  the  immunogenicity  of  an  -14  kDa  protein 
was  low  relative  to  the  immunogenicity  of  an  -38  kDa  protein,  a 
statistically  significant  difference  between  Leishmania- infected 
donors  and  controls  was  found  (P  <  0.05,  Student’s  t  test). 
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Figure  2.  Immunogenicity  to  P.  papcitcisi  salivary  gland  proteins  among  plasma  donors.  (A)  Immunogenicity  of  all  P.  papatasi  salivary  gland 
proteins  by  Egyptian  blood  plasma  donor  antibody  specificities  from  three  regions  (CA,  AW,  and  NS)  using  consensus  estimates  of  MW  on 
immunoblot  images  from  three  independent  observers.  (B)  Immunogenicity  of  all  P.  papatasi  salivary  gland  proteins  by  Jordanian  blood  plasma 
donor  antibody  specificities  from  three  regions  (AM,  MA,  and  SW).  Ecological  conditions  in  the  regional  key  are  the  same  as  those  conditions 
outlined  in  Table  1.  (C)  Immunogenicity  of  all  P.  papatasi  salivary  gland  proteins  of  Iraq-based  US  personnel  and  blood  plasma  donor  antibody 
specificities  from  individuals  differentially  exposed  to  sand  fly  bites  and  Leishmania  parasites.  Exposure  criteria  are  the  same  as  those  criteria  outlined 
in  Table  1.  Immunogenicity  is  calculated  as  (probability  of  donor  plasma  antibody  reaction  of  protein  among  a  regional  group)  x  (average  of  band 
intensity  among  that  regional  group)  =  a.u.  *Significant  difference(s)  in  host  plasma  antibody  specificity  for  protein(s)  at  these  MWs  between  sand  fly- 
endemic  and  control  regional  sites  or  between  sand  fly  bite-exposed  military  personnel  and  unexposed  control  donors  ( P  <  0.05). 


Despite  significant  relative  differences  of  raw  immunoblot 
scores  to  the  control  region  for  Egyptian  donors,  immunoge¬ 
nicity  values  for  salivary  proteins  at  ~26  and  44  kDa  were  low 
compared  with  those  values  of  —12,  30,  and  38  kDa  proteins. 
Both  the  26  and  44  kDa  proteins  were  found  in  close  relative 
proximity  to  one  or  more  proteins  of  similar  MW  (Figure  IB). 
Because  each  donor  plasma  profile  displayed  differential  spec¬ 
ificity  for  majorly  expressed  salivary  proteins,  manual  scoring 
may  have  been  imprecise  for  proteins  at  those  two  MWs. 
Therefore,  we  took  additional  steps  to  more  thoroughly  char¬ 


acterize  both  the  majorly  expressed  salivary  proteins  and  the 
global  single  human  donor  plasma  antibody  specificity  profiles. 

Sand  fly  salivary  protein  identification.  To  more  thoroughly 
differentiate  and  characterize  abundantly  expressed  P.  papatasi 
salivary  proteins,  10  bands  were  excised  from  SGS  segregated 
by  denaturing  PAGE  (~12,  14,  26,  28,  30,  38,  40,  42,  44,  and 
46  kDa).  MS  analyses  were  performed  on  the  purified  proteins 
after  trypsin  digestion.  Resultant  spectra  were  compared 
against  all  available  NCBI  protein  database  sequences  for 
Phlebotomus  and  Lutzomyia  spp.33  Table  3  outlines  the  top 
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Table  3 


Identification  of  majorly  expressed  P.  papatasi  salivary  proteins  by  MS 


P.  papatasi  salivary 
protein  (estimated 
MW)*/rank 

Unused 

Sequence 
coverage t  (%) 

Accession 

number 

Protein  name 

Species 

No.  of  unique 
peptides 

Confidence  of 
identification^ 

Local  FDR  for 
protein  (%) 

12  kDa 

1 

50.41 

80 

gi  15963505 

12  kDa  salivary  protein  precursor 

P.  papatasi 

56 

~>  le-50 

7e-5 

2 

25 

57.8 

gi  15963507 

14  kDa  salivary  protein  precursor 

P.  papatasi 

16 

~>  le-25 

7.3e-5 

3 

6 

27.2 

gi 15963515 

32  kDa  salivary  protein  precursor 

P.  papatasi 

4 

~>  le-6 

7.54e-5 

4 

5.56 

55.3 

gi 112497698 

14.2  kDa  salivary  protein 

P.  duboscqi 

6 

~>  le-5 

7.66e-5 

5 

4 

72.1 

gi 112497317 

13.7  kDa  salivary  protein 

P.  duboscqi 

22 

~>  le-4 

7.85e-5 

6 

2 

30.5 

gi 112497287 

14.2  kDa  salivary  protein 

P.  duboscqi 

1 

0.01 

8.04e-5 

14  kDa 

1 

33.3 

86.1 

gi  335954166 

SP15  protein 

P.  papatasi 

50 

~>  le-33 

8.24e-5 

2 

11.18 

86.1 

gi  335954098 

SP15  protein 

P.  papatasi 

44 

~>  le-11 

8.57e-5 

3 

6 

28.2 

gi  15963507 

14  kDa  salivary  protein  precursor 

P.  papatasi 

3 

~>  le-6 

8.66e-5 

4 

4.05 

90.9 

gi 112497496 

14.5  kDa  salivary  protein 

P.  duboscqi 

39 

~>  le-36 

8.76e-5 

5 

4 

10.6 

gi  15963515 

32  kDa  salivary  protein  precursor 

P.  papatasi 

2 

~>  le-4 

8.86e-5 

6 

2 

39.4 

gi  112361953 

14.4  kDa  salivary  protein 

P.  duboscqi 

3 

0.01 

9.05e-5 

7 

2 

32.1 

gi 157361529 

40S  ribosomal  protein  S12-like  protein 

P.  papatasi 

1 

0.01 

9.05e-5 

8 

0.26 

86.1 

gi  335954106 

SP15  protein 

P.  papatasi 

40 

0.55 

9.55e-5 

26  kDa 

1 

30.74 

71.3 

gi  15963511 

28  kDa  salivary  protein  precursor 

P.  papatasi 

36 

~>  le-30 

8.78e-5 

2 

28.47 

70 

gi  15963513 

30  kDa  salivary  protein  precursor 

P.  papatasi 

40 

~>  le-28 

8.85e-5 

28  kDa 

1 

30.78 

57.5 

gi  15963511 

28  kDa  salivary  protein  precursor 

P.  papatasi 

26 

~>  le-30 

4.69e-5 

2 

23.26 

60.5 

gi 15963513 

30  kDa  salivary  protein  precursor 

P.  papatasi 

18 

~>  le-23 

4.72e-5 

30  kDa 

1 

28.72 

56.3 

gi  76589378 

29  kDa  salivary  protein  precursor 

P.  papatasi 

29 

~>  le-28 

9.06e-5 

2 

12 

37.8 

gi 15963515 

32  kDa  salivary  protein  precursor 

P.  papatasi 

8 

~>  le-12 

9.25e-5 

38  kDa 

1 

113.02 

92 

gi  10443907 

Salivary  apyrase 

P.  papatasi 

114 

~>  le-113 

4.94e-5 

2 

6.06 

74.1 

gi  112496715 

35.8  kDa  salivary  protein 

P.  duboscqi 

27 

~>  le-6 

5.21e-5 

3 

2.73 

32 

gi 299829440 

32.9  kDa  salivary  protein 

P.  sergenti 

7 

~>  le-3 

5.31e-5 

40  kDa 

1 

72.21 

80.8 

gi  15963517 

42  kDa  salivary  protein  precursor 

P.  papatasi 

70 

~>  le-72 

5.08e-5 

2 

9.68 

51.1 

gi 112497236 

Yellow-related  protein 

P.  duboscqi 

29 

~>  le-9 

5.15e-5 

3 

7.52 

69.1 

gi  10443907 

Salivary  apyrase 

P.  papatasi 

6 

~>  le-72 

5.24e-5 

42  kDa 

1 

120.09 

89.6 

gi  15963517 

42  kDa  salivary  protein  precursor 

P.  papatasi 

110 

~>  le-120 

1.04e-3 

2 

21.3 

73 

gi  15963519 

44  kDa  salivary  protein  precursor 

P.  papatasi 

13 

~>  le-21 

1.09e-3 

3 

15.46 

65.6 

gi 112497236 

Yellow-related  protein 

P.  duboscqi 

46 

~>  le-15 

l.lle-3 

4 

4 

55.2 

gi 299829410 

Yellow-related  salivary  protein 

P.  duboscqi 

19 

~>  le-4 

1.12e-3 

44  kDa 

1 

91.06 

85.8 

gi  15963517 

42  kDa  salivary  protein  precursor 

P.  papatasi 

100 

~>  le-90 

6.72e-7 

2 

81.47 

80.3 

gi 15963519 

44  kDa  salivary  protein  precursor 

P.  papatasi 

97 

~>  le-81 

6.72e-7 

3 

15.45 

61.8 

gi  112497236 

Yellow-related  protein 

P.  duboscqi 

49 

~>  le-15 

6.75e-7 

4 

6.13 

59.7 

gi 112497202 

Yellow-related  protein 

P.  duboscqi 

41 

~>  le-6 

6.76e-7 

46  kDa 

1 

130.07 

98.3 

gi  15963519 

44  kDa  salivary  protein  precursor 

P.  papatasi 

188 

~>  le-130 

l.le-5 

2 

43.93 

77.2 

gi 15963517 

42  kDa  salivary  protein  precursor 

P.  papatasi 

41 

~>  le-43 

l.le-5 

3 

8.89 

74.7 

gi 112497202 

Yellow-related  protein 

P.  duboscqi 

72 

~>  le-8 

l.lle-5 

4 

6.02 

21.9 

gi  76446619 

Yellow-related  salivary  protein 

P.  perniciosus 

3 

~>  le-6 

1.12e-5 

5 

4.29 

35.1 

gi  10443907 

Salivary  apyrase 

P.  papatasi 

3 

~>  le-4 

1.12e-5 

6 

2.26 

50.1 

gi 112497236 

Yellow-related  protein 

P.  duboscqi 

26 

~>  le-3 

1.13e-5 

7 

2 

38.7 

gi  52001011 

Yellow-related  salivary  protein 

P.  papatasi 

1 

0.01 

1.14e-5 

*  Based  on  denatured  protein  bands  excised  for  MS  analysis  (Figure  IB). 

t  Percent  representation  of  a  protein  based  on  unique  peptide  spectral  assignments  to  that  protein. 

t  Negative  loglO  of  the  unused  ProtScore,  which  is  calculated  using  only  peptides  from  spectra  that  have  not  been  used  by  other  proteins;  a  lower  value  indicates  greater  confidence  of  protein 
identification-based  additive  assignments  of  unique  peptides. 


protein  matches,  ordered  by  rank  of  best  local  match  to  our 
database,  for  spectra  of  each  protein  band.  All  top-ranked 
matches  were  against  published  protein  sequences  from  the 
P.  papatasi  Israeli  colony  strain.  Despite  seemingly  clear  and 
unique  in-gel  MW  differences  for  proteins  in  the  26  and  28  kDa 
bands  as  well  as  those  proteins  in  the  40,  42,  44,  and  46  kDa 
bands,  multiple  high-confidence  protein  matches  were  found 
for  each  analyzed  band.  Moreover,  given  the  high  homology 
in  primary  sequence  for  PpSP28  and  PpSP30  (top  matches  for 
the  ~26  and  28  kDa  bands)  as  well  as  between  PpSP42  and 
PpSP44  (top  matches  for  the  ~40,  42,  44,  and  46  kDa  bands), 
identification  was  based  predominantly  on  spectral  matches  to 


domains  exclusively  unique  to  each  protein.  Within  the  range 
of  -26-30  kDa,  we  distinguished  three  primary  salivary  pro¬ 
teins  that  corresponded  to  differential  donor  antibody  specific¬ 
ities  (Figure  IB),  where  MS  matches  to  the  -26  and  28  kDa 
proteins  were  equally  confident  between  PpSP28  and  PpSP30. 
The  band  at  -30  kDa  matched  best  to  a  29  kDa  salivary  protein 
precursor  (gi|76589378). 

Clustering  of  host  antibody  anti-salivary  protein  profiles. 

K-means  clustering  was  performed  on  raw  immunoblot  score 
data,  and  therefore,  the  complete  global  antibody  specificity 
profile  for  each  individual  donor  was  simultaneously  considered. 
This  approach  helped  to  better  characterize  the  contribution 


930 


GERACI  AND  OTHERS 


Table  4 

Principle  salivary  proteins  in  K-means  clusters  (estimated  MW) 


12  kDa 

26  kDa 

30  kDa 

38  kDa 

44  kDa 

MENA  clusters 

A 

B 

+ 

- 

+ 

+ 

- 

C 

+ 

+ 

+ 

+ 

- 

D 

+ 

+ 

+ 

+ 

+ 

US  military  clusters 

12  kDa 

14  kDa 

26  kDa 

30  kDa 

38  kDa 

E 

- 

- 

- 

- 

- 

F 

- 

- 

- 

- 

+ 

G 

- 

+ 

+ 

+ 

+ 

H 

+ 

+ 

- 

+ 

- 

Cluster  characteristics  are  defined  as  antibody  recognition  at  the  listed  protein  MW  for 
>  50%  of  donors.  Results  details  descriptions  of  clusters. 


that  regional  and  ecological  characteristics  played  on  host- 
acquired  immunity  to  sand  fly  salivary  proteins  by 
distinguishing  common  patterns  within  the  large  complex 
dataset.  K-means  clustering  was  performed  first  with  Egyptian 
and  Jordanian  donor  immunoblot  data  using  k  =  6  (one  cluster 
for  each  significant  protein  plus  a  null  control)  and  again  with 
k  =  5  (the  number  of  significant  proteins).  No  distinguishing  or 
consistently  similar  data  patterns  were  observed  using  these 
numbers  of  clusters  (data  not  shown).  However,  on  using  k  = 
4  clusters  (Supplemental  Figure  3),  data  patterns  emerged  that 
were  both  consistent  among  donors  within  each  single  cluster 
and  significantly  divergent  between  clusters  by  cumulative 
immunoblot  intensity  scores  for  the  five  immunogenic  salivary 
proteins  (P  <  0.05,  Student’s  t  test).  Similarly,  in  clustering 
immunoblot  data  from  the  US  military  donors,  the  use  of  k  = 
4  clusters  generated  cluster  profiles  with  distinct  patterns  of 
specificity  against  the  -38,  30,  26,  14,  and  12  kDa  salivary 
proteins  (Supplemental  Figure  4).  Each  cluster  was  also  signif¬ 
icantly  divergent  in  this  dataset  by  cumulative  immunoblot 
intensity  scores  for  this  different  set  of  five  distinguishing  sali¬ 
vary  proteins  (P  <  0.05,  Student’s  t  test).  Table  4  outlines  the 
characteristic  plasma  antibody  specificity  patterns  for  each  of 
the  five  immunogenic  salivary  protein  sets. 

Clusters  A-D  pertain  to  Egyptian  and  Jordanian  resident 
datasets.  Cluster  A  is  a  null  response  cluster  (populated  pre¬ 
dominantly  by  control  region  donors).  Cluster  B  is  a  cluster 
with  antibodies  specific  to  -12,  30,  and  38  kDa  salivary  pro¬ 
teins  (minor  difference  between  the  numbers  of  donors  from 
each  region).  Cluster  C  is  a  cluster  with  predominant  speci¬ 
ficity  to  -12,  30,  38  and  26  kDa  proteins  (again,  no  substan¬ 
tial  differences  in  regional  donor  numbers).  Cluster  D  is  a 
cluster  with  specificity  to  -12,  26,  30,  38,  and  44  kDa  salivary 
proteins  (substantially  more  donors  residing  in  sand  fly-  and 
leishmaniasis-endemic  regions  compared  with  controls).  Only 
the  percentages  of  Egyptian  CA  and  AW  donors  displayed  a 
significant  inverse  Pearson  correlation  across  clusters  (P  < 
0.05,  r  =  -0.94).  However,  the  percentages  of  donors  from 
both  AW  and  NS  trended  positively  across  clusters,  whereas 
the  percentages  of  CA  donors  trended  negatively  (Figure  3A). 
For  Jordanian  donors,  a  significant  positive  correlation  was 
found  between  AM  and  MA  donors  across  clusters  (P  <  0.05, 
r  =  0.96),  and  no  correlations  were  observed  between  AM 
and  SW  (Figure  3B).  No  correlations  existed  across  clusters 
between  all  donor  immunoreactivity  profiles  combined  by 
ecotope  (Figure  3C). 

Clusters  E-H  pertain  to  the  US  military  donor  dataset. 
Cluster  E  is  a  null  response  cluster  (populated  predominantly 
by  control  region  donors).  Cluster  F  is  a  cluster  with  antibod¬ 


ies  specific  almost  exclusively  to  the  -38  kDa  salivary  protein 
(a  substantial  difference  in  representation  from  each  exposure 
group,  where  the  greatest  percentage  of  individuals  was  from 
the  Leishmania-'mte.c\.e.&  group).  Cluster  G  is  a  cluster  with 
predominant  specificity  to  -14,  30,  26,  and  38  kDa  proteins 
(again,  a  substantial  difference  in  representation  from  each 
exposure  group,  where  the  greatest  percentage  was  from  the 
bitten  but  not  infected  group).  Cluster  H  is  a  cluster  with 
specificity  to  -12,  14,  and  30  kDa  salivary  proteins  (a  cluster 
populated  largely  of  control  donors).  No  statistically  signifi¬ 
cant  correlations  existed  between  exposure  groups  across 
clusters  E-H.  The  distribution  of  donors  favored  a  single  clus¬ 
ter,  where  one  anti-salivary  gland  protein  specificity  profile 
was  represented  predominantly  by  donors  of  a  single  expo¬ 
sure  type  group  (Figure  3D). 

Anti-SGS  antibody  isotype  assessments.  To  further  under¬ 
stand  the  acquired  antibody-based  immunity  to  sand  fly  sali¬ 
vary  proteins,  the  six  most  representative  donor  plasmas  from 
each  of  the  MENA  regional  resident  study  sites  (associated 
with  each  of  the  four  k-means  clusters)  were  assessed  for 
specific  anti-SGS  Ig  isotype  classifications  (Figure  4  and  Sup¬ 
plemental  Figure  5).  Compared  with  concentrations  for  all 
other  antibody  isotypes,  significantly  greater  cumulative  con¬ 
centrations  of  anti-SGS  IgG4  were  present  in  donor  plasma 
(P  <  0.05)  (Figure  4).  Notably,  when  data  for  all  Ig  isotypes 
was  analyzed  cumulatively,  a  significant,  albeit  weak,  positive 
Pearson  correlation  was  observed  across  the  four  clusters 
(P  =  0.01,  r  =  0.23).  When  each  isotype  was  analyzed  sepa¬ 
rately,  IgGl,  -2,  and  -4  and  IgE  concentrations  were  positively 
correlated  across  the  clusters  (P  =  0.01-0.05,  r  =  0.42-0.54). 
Only  IgG3  did  not  correlate  with  clusters,  likely  as  a  result  of 
such  low  relative  concentrations.  These  correlations  logically 
follow  the  defining  characteristics  of  each  cluster’s  immuno¬ 
reactivity  profile,  wherein  progressively,  there  is  greater  spec¬ 
ificity  to  more  salivary  proteins  in  clusters  A-D.  Evidence 
exists  for  a  correlation  between  anti-mosquito  salivary  IgG4 
and  IgE  expression,  which  contributes  to  host  type  1  hyper¬ 
sensitivity.40  Here,  anti-/',  papatasi  SGS  IgG4  and  IgE  concen¬ 
trations  did  correlate  modestly  (P  =  0.02,  r  =  0.5).  Whether 
this  correlation  is  actually  related  to  the  establishment  of  a 
true  hypersensitivity  among  donors  residing  in  our  ecotopes 
of  study  remains  unknown. 

Cellular  proliferation  analysis.  To  assess  donor  cellular 
responsiveness  to  P.  papatasi  salivary  proteins,  peripheral  blood 
mononuclear  cells  from  a  random  subset  of  donors  from  each 
regional  location  and  exposure  group  were  assessed  in  a  prolif¬ 
eration  assay.  Monocytic  cells  were  adherence-selected  from 
cryopreserved  donor  leukocytes  and  differentiated  to  dendritic 
cells  by  cytokine  stimulation.  Donor  MDDCs  were  exposed  to 
SGS  1  day  before  coculture  with  remaining  autogenic  leuko¬ 
cytes.  Cellular  proliferation  between  SGS  pre-stimulated  and 
unstimulated  cocultures  was  measured  by  incorporation  of  tri- 
tiated  thymidine  (Figure  5).  MDDCs  of  donors  residing  in  sand 
fly-endemic  regions  elicited,  on  average,  less  proliferation  after 
SGS  exposure  than  those  MDDCs  from  control  region  donors 
(Figure  5A).  However,  only  the  stimulation  ratios  compared 
between  AM  and  MA  site  donors  were  significantly  different 
(Student’s  t  test,  P  <  0.05).  Moreover,  Iraq-deployed  but  unin¬ 
fected  US  military  donors  exhibited  less  SGS-induced  prolifer¬ 
ation  compared  with  control  donors  (Figure  5B). 

Cytokine  secretion  profiles.  To  better  understand  cellular 
sensitization  of  immunologically  naive  individuals  to  P.  papatasi 
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Figure  3.  K-means  clustering  of  total  host  plasma  antibody  specificity  profiles  to  P.  papatasi  whole-salivary  gland  protein  sonicates.  Scores  for 
human  host  plasma  specificity  to  sand  fly  salivary  proteins  were  separated  by  national  and  ecological  region  groups  and  then  clustered  by  the  k-means 
method  across  1,000  iterations  using  Cluster  v3.0,  with  k  =  4  clusters  and  the  Euclidean  distance  metric.  (A)  Total  percentages  of  blood  plasma  donors 
from  each  Egyptian  regional  site  assigned  to  each  anti-SGS  antibody  specificity  profile  cluster  by  k-means  clustering.  (B)  Total  percentages  of  blood 
plasma  donors  from  each  Jordanian  regional  site  assigned  to  each  anti-SGS  antibody  specificity  profile  cluster  by  k-means  clustering.  (C)  Total 
percentages  of  blood  plasma  donors  from  combined  regional  sites  (grouped  by  defining  ecological  characteristics)  (Table  1)  assigned  to  each  anti- 
SGS  antibody  specificity  profile  cluster  by  k-means  clustering.  (D)  Total  percentages  of  blood  plasma  donors  from  US  military  personnel  differentially 
exposed  to  sand  fly  bites  and  Leishmania  parasites  assigned  to  each  anti-SGS  antibody  specificity  profile  cluster  by  k-means  clustering.  *K-means 
cluster  profile  characteristics  are  outlined  in  Table  4.  t  A  significant  Pearson  correlation  was  found  between  international  donors  from  these  sites. 


salivary  proteins,  proliferation  culture  supernatants  of  Iraq- 
based  military  donors  exposed  to  SGS  were  analyzed  to  quan¬ 
tify  21  secreted  cytokine  concentrations.  Interleukin-8  (IL-8), 
IL-17,  and  interferon-y-induced  protein-10  (IP-10)  of  SGS- 
stimulated  donor  cultures  were  significantly  greater  in  the 
supernatants  of  control  exposure  group  donors  than  the  super¬ 
natants  of  cells  from  the  Leishmania-mfcctcd  donors  (P  <  0.05, 
Student’s  t  test)  (Figure  6). 


DISCUSSION 

We  present  here  an  examination  of  human  host  responses 
to  P.  papatasi  salivary  proteins.  The  overarching  goal  of  this 
study  was  to  understand  the  role  of  natural  exposure  to  wild 
sand  flies  by  characterizing  host-acquired  immunity  dynamics 
with  analyses  of  circulating  antibody  specificities  and  isotypes 
to  whole-sand  fly  SGSs  as  well  as  the  memory  response  of 
among  leukocytes  from  those  same  donors.  Together,  these 


data  serve  to  further  our  understanding  of  host  responses  to 
sand  fly  vector  salivary  proteins.  Sites  were  selected  in  both 
Egypt  and  Jordan  by  presumed  differential  exposure  of  resi¬ 
dents  to  wild  P.  papatasi  populations.  In  the  context  of  human 
antibody  specificities  against  SGS  proteins,  a  starker  contrast 
between  ecotope  sites  was  observed  among  Egyptian  sites 
compared  with  sites  in  Jordan  (Figure  2A  and  B  and  Table  2). 
All  collection  sites  in  Egypt  were  located  farther  apart  from  one 
another  compared  with  Jordanian  collection  sites  (100+  miles 
versus  ~20  miles  between  sites,  respectively),  with  far  less 
intersite  human  traffic.  Such  proximity  and  easier  travel 
throughout  Jordan  may  be  a  factor  contributing  to  SGS  reac¬ 
tivity  of  donors  residing  in  control  regions.  However,  in  addi¬ 
tion  to  potential  donor  travel  to  endemic  P.  papatasi  regions, 
the  possibility  of  exposure  to  non-sand  fly  hematophagous 
insects,  generating  antibodies  with  specificity  to  conserved  pro¬ 
tein  domains,  should  not  be  discounted.41 

Our  data  for  MENA  residents  compare  well  with  the  data  in 
the  work  by  Marzouki  and  others,  which  analyzed  plasma 
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Figure  4.  Isotyping  of  anti-SGS  human  host  plasma  antibodies.  For  each  k-means  cluster,  six  plasmas  most  representative  of  that  k-means 
cluster  immunoreactivity  profile  and  plasmas  from  individuals  residing  in  each  of  the  Middle  Eastern  regions  were  isotyped  for  total  anti-SGS 
antibodies  using  custom  indirect  ELISA.  Individual  donor  results  for  each  antibody  isotype  may  be  found  in  Supplemental  Figure  5.  ^Clusters  are 
like  those  clusters  described  in  Table  4.  tAnti-SGS  antibody  concentrations  were  assessed  by  secondary  antibody  detection  of  a  known  concen¬ 
tration  of  purified  Ig  isotype  standards,  plated  in  dilution  series  as  primary  antigen,  and  fitted  to  the  resultant  standard  curve;  final  values  were 
averaged  per  cluster  across  all  regional  donors. 


antibody  specificity  of  Tunisian  donors  to  P.  papatasi  salivary 
proteins.  Marzouki  and  others12  found  that  the  majority  of 
their  donors  exhibited  positive  reactivity  to  salivary  proteins  at 
~12, 15,  21,  28,  30,  36,  and  44  kDa.  Sand  fly  vectors  in  Tunisia 
are  phylogenetically  distinct  from  sand  flies  endemic  to  Egypt, 
Jordan,  and  Israel,  which  are  part  of  a  separate,  albeit  some¬ 
what  heterogeneous,  population.42  The  laboratory-reared 
P.  papatasi  used  in  the  present  study  originated  from  Israel. 

Glycoinositolphospholipids  (GIPLs)  on  the  surface  of 
Trypanosoma  spp.  have  been  shown  to  alter  the  salivary  pro¬ 


file  of  their  insect  vector,43  invoking  the  possibility  that  cross¬ 
strain  variability  among  Leishmania  spp.  might  affect  vector 
saliva  profiles  and  consequently,  change  human  host  anti- 
salivary  protein  antibody  reactivity.  Failure  of  Egyptian 
and  Jordanian  donor  plasma  antibodies  to  recognize  the  sali¬ 
vary  protein  at  ~14  kDa  (Figure  1),  which  was  identified  by 
Marzouki  and  others,12  may  be  because  of  the  differential 
expression  of  this  protein  in  distinct  sand  fly  populations. 
However,  US  military  personnel  exposed  to  sand  fly  bites  in 
Iraq  displayed  specificity  to  the  protein  at  ~14  kDa.  The  fact 
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Figure  5.  Mixed  leukocyte  proliferation  in  vitro  of  peripheral  blood  donors  after  stimulation  of  MDDCs  with  P.  papatasi  SGS.  (A)  Represen¬ 
tative  sampling  of  regional  Egyptian  and  Jordanian  resident  donor  cellular  proliferations  when  exposed  to  SGS-stimulated  autogeneic  MDDCs. 
Proliferation  was  measured  using  coculture  with  tritiated  thymidine  for  24  hours  and  (I-particle  decay  counts  per  minute.  Values  presented  are 
ratios  of  stimulated  over  unstimulated  assays  normalized  to  a  single  positive  control  donor  proliferation  run  concurrently.  (B)  Proliferation  data 
for  US  military  personnel  differentially  exposed  to  sand  fly  bites  and  Leishmania  parasites.  ““Significant  differences  by  Student's  t  test  (P  <  0.05). 
#Stimulation  index  is  the  ratio  of  (I-counts  (3H-thymidine  uptake)  between  unstimulated  and  SGS  pre-stimulated  MDDCs  and  mixed  leuko¬ 
cyte  cocultures. 
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Figure  6.  Secreted  cytokines/chemokines  of  SGS-stimulated 
MDDCs.  Culture  supernatant  expression  of  (A)  IL-8,  (B)  IL-17,  and 
(C)  IP-10  secreted  by  in  vitro  MDDCs  from  differential  sand  fly-  and 
Leishmania-ex posed  US  military  personnel  blood  donors  after  stimu¬ 
lation  with  SGS.  *Data  are  MFI  +  SD  after  normalization  of  SGS- 
stimulated  MDDC  supernatants  to  paired  unstimulated  supernatant 
values.  **Significant  differences  by  Student’s  t  test  (P  <  0.05). 


that  this  previously  naive  group  produced  antibodies  specific 
to  this  protein  and  endemic  residents  did  not  highlights  the 
difference  between  chronic  exposure  to  P.  papatasi  bites  ver¬ 
sus  initial  contact  with  new  antigens  in  immunologically  naive 
persons.  Indeed,  the  study  by  Marzouki  and  others12  used 
blood  samples  collected  from  children  between  the  ages  of  6 
and  12  years.  The  young  age  of  donors  in  that  study,  likely 
lacking  the  prolonged  history  of  sand  fly  bite  exposure 
encountered  by  adult  endemic  residents  in  our  current  work, 
may  account  for  the  differential  specificities  observed  against 
the  ~14  kDa  salivary  protein.  The  fact  that  the  immunogenic- 
ity  to  this  protein  among  US  personnel  was  significantly  less 
within  the  Leishmania- infected  donor  group  compared  with 
the  sand  fly  bite  only  donors  and  controls  suggests  a  possible 
role  of  the  parasite  in  diminishing  host  antibody  production 
against  the  ~14  kDa  protein  or  a  possible  impact  on  sand  fly 
expression  or  post-translational  modification  of  the  protein 
before  host  contact. 

The  lack  of  expression  of  an  ~21  kDa  protein  (Figure  IB) 
by  our  Israeli  colony  of  P.  papatasi  could  be  because  of  differ¬ 
ential  salivary  gene  expression  between  variant  vector  popu¬ 
lations,  especially  given  the  seasonal  and  regional  differences 
in  expression  of  select  salivary  protein  transcripts.15  Moreover, 
differences  in  post-translational  modifications  of  salivary  pro¬ 
teins,  such  as  glycosylation  or  phosphorylation  patterns,44 
between  the  Tunisian  and  Israeli  colonies  may  have  led  to 
differences  in  protein  identification  and  even  antibody  specific¬ 
ities.  Both  studies  used  manual  scoring  of  immunoblots,  and 
designations  of  positive  specificities  at  particular  MWs  may  not 
be  consistent  between  studies  (i.e.,  what  our  observers  agreed 
was  specificity  at  26  and  38  kDa  may  have  been,  for  Marzouki 
and  others12  observers,  specificity  at  28  and  36  kDa).  It  should 
also  be  noted  that  Marzouki  and  others12  identify  their  linear¬ 
ized  21  kDA  salivary  protein  as  PpSP42;  we  also  identified  this 
protein,  corresponding  to  linearized  salivary  proteins  in  the 
-4CM16  kDa  range. 

To  assess  the  possible  biological  significance  of  immuno¬ 
genic  profile  characteristics,  P.  papatasi  salivary  proteins  were 
identified  by  MS  analysis.  Proteins  determined  to  be  ~26  and 
28  kDa  by  denaturing  PAGE  matched  with  approximately 
equivalent  coverage  and  confidence  to  PpSP28  and  PpSP30. 
The  mature  MWs  for  PpSP28  (gi|  15963511)  and  PpSP30 
(gi|  15963513)  are  27.2  and  27.8  kDa,  respectively,  and  both 
proteins  share  82.1%  primary  sequence  identity  by  pairwise 
alignments,  with  multiple  homologous  domains.  This  infor¬ 
mation  suggests  that  immunoblot  scores  for  proteins  in  the 
26-28  kDa  range  may  be  for  one  or  both  proteins  by  cross- 
immunoreactivity.  PpSP28  and  PpSP30  transcripts  were  dif¬ 
ferentially  expressed  among  wild-caught  P.  papatasi  across 
sites  and  time.15  Hence,  exposure  to  sand  fly  bites  at  vari¬ 
ous  seasons  throughout  a  year  may  influence  global  immune 
responsiveness  to  P.  papatasi.  Although  all  blood  samples  from 
the  MENA  region  in  this  study  were  collected  during  a  3-week 
window  of  late  summer  in  2007,  circulating  anti-salivary  anti¬ 
bodies  may  be  a  result  of  prior  exposures.  The  functional  pre¬ 
dictions  for  PpSP28  and  PpSP30  implicate  their  role  in  binding 
thromboxane  A,41  a  host  factor  involved  in  vasoconstriction 
and  prothromobosis,  a  process  that  would  be  advantageous  to 
stifle  the  promotion  of  effective  vector  blood  feeding. 

Our  results  revealed  that  one  majorly  expressed  salivary 
protein  of  ~30  kDa  is  immunodominant  to  other  salivary 
proteins  among  MENA  donors  (Figure  2  and  Table  2). 
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Another  report  implicated  a  ~30  kDa  salivary  protein  as 
strongly  immunodominant,13  and  Marzouki  and  others11 
recently  postulated  that  the  strongest  degree  of  antibody  spec¬ 
ificity  among  pre-exposed  human  and  murine  hosts  was  against 
a  32  kDa  P.  papatasi  salivary  protein  (PpSP32,  gi|  15963515). 
We  analyzed  an  SGS  protein  band  of  approximately  the  same 
MW  by  MS  (Figure  IB).  The  best  match  from  among  all  pub¬ 
lished  phlebotomine  sand  fly  protein  sequences33  was  PpSP29 
(gi|  76589378)  (Table  2).  PpSP32  emerged  as  a  secondary 
match  for  the  same  protein  gel  band.  PpSP29  and  PpSP32 
share  little  homology  between  primary  sequences  by  pair¬ 
wise  alignment  (data  not  shown).  Most  P.  papatasi  protein 
sequences  in  our  local  database  were  those  sequences  pub¬ 
lished  from  work  with  the  Israeli  colony,  the  same  colony 
used  throughout  our  study.  Marzouki  and  others12  performed 
experiments  using  a  Tunisian  strain  of  P.  papatasi  but  per¬ 
formed  identifications  using  the  same  approximate  database 
of  protein  sequences.  It  follows  logically  that  our  MS  analysis 
yielded  matches  with  higher  degrees  of  confidence  compared 
with  the  study  by  Marzouki  and  others,12  which  used  a  pop¬ 
ulation  of  differential  phylogenetic  origin42  from  the  com¬ 
mon  protein  database.  Abdeladhim  and  others41  recently 
performed  an  exhaustive  analysis  of  the  Tunisian  strain 
P.  papatasi  salivary  proteome  through  translated  Basic  Local 
Alignment  Search  Tool  (BLAST)  alignments  of  RNA  tran¬ 
scripts  to  the  same  common  protein  sequence  database. 
PpSP29  was  found  to  display  homology  to  insect  allergen 
5  (AG5),  a  protein  expressed  in  most  hematophagous  insects. 
Immunoreactivity  to  PpSP29,  given  its  membership  in  the  pan¬ 
vector  AG5  protein  family,  could  be  generated  by  exposure  to 
other  insect  vectors  because  of  conserved  immunogenic  epi¬ 
topes.  PpSP32  was  predicted  to  be  a  potential  mucin  based  on 
similarities  of  glycosylation  sites  and  thus  far,  has  been  shown 
to  be  exclusively  produced  by  sand  flies.45  Exposure  to  mucin¬ 
like  PpSP32  might  stimulate  T-cell  proliferation  and  cytokine 
production,  which  has  been  known  to  occur  in  response  to 
other  true  mucins.46  Given  the  agreed  importance  of  proteins 
at  the  30-kDa  range,  future  analysis  using  recombinant 
PpSP29  and  PpSP32  will  be  beneficial  to  ascertain  which  pro¬ 
tein  is  the  true  immunodominant  salivary  antigen. 

Conversely,  for  Iraq-based  US  military  donors,  the  most 
immunodominant  P.  papatasi  salivary  protein  was  one  of 
~38  kDa  (Figure  2  and  Table  2).  Marzouki  and  others1112  also 
noted  that  Tunisian  donor  antibodies  displayed  high  specific¬ 
ity  to  this  protein  (referred  to  as  an  ~36  kDa  protein)  but  not 
the  same  degree  as  the  -30/32  kDa  protein(s).  The  -38  kDa 
protein  matched  best  to  a  salivary  apyrase  (gi|  10443907). 
Apyrases  hydrolyze  both  adenosine  triphosphate  (ATP)  and 
adenosine  diphosphate  (ADP)  to  adenosine  monophosphate 
(AMP),  thus  destroying  an  important  physiologic  stimulus  of 
host  platelet  aggregation  released  from  damaged  tissues  and 
blood  cells.  Apyrases  have  been  found  in  most  every  hema¬ 
tophagous  arthropod.47  The  stark  dichotomy  between  two 
donor  groups  (presumably  chronically  exposed  regional  resi¬ 
dents  and  previously  immunologically  naive  US  military)  in 
terms  of  reactivity  to  the  -30  and  38  kDa  proteins  may 
suggest  a  population  scale-acquired  immune  specificity  for  a 
protein  that  could  be  more  specific  to  the  Phelobotomus 
genus  (-30  kDa),  versus  one  that  is  more  encountered  from 
insect  bites  occurring  in  the  United  States  (-38  kDa),  or  at 
least  to  which  immunologically  naive  individuals  more  readily 
produce  antibodies  against  after  initial  exposures. 


We  used  k-means  clustering  to  identify  patterns  of  plasma 
antibody  immunoreactivity  to  SGS  proteins  across  donors. 
Antibody  specificity  to  the  five  salivary  proteins  that  elicited 
significantly  different  immunoreactivity  among  MENA- 
endemic  donors  compared  with  control  region  data  was  the 
principle  factor  distinguishing  clusters  A-D  (Table  4  and  Sup¬ 
plemental  Figure  3).  Among  Egyptian  donors,  CA  and  AW 
displayed  a  significant  inverse  Pearson  correlation  across  clus¬ 
ters  (P  <  0.05,  r  =  -0.94).  No  significant  inverse  correlation 
was  observed  between  CA  and  NS,  despite  an  inverse  trend 
(Figure  3A).  For  Jordanian  donors,  a  significant  positive  cor¬ 
relation  was  found  between  AM  and  MA  donors  across  clus¬ 
ters  ( P  <  0.05,  r  =  0.96),  and  no  correlations  were  observed 
between  AM  and  SW  (Figure  3B).  The  shorter  intersite  dis¬ 
tances  and  increased  travel  of  residents  in  Jordan  may  have 
led  to  more  homogeneous  exposure  conditions  compared 
with  the  more  distant  and  isolated  Egyptian  sites,  accounting 
for  the  negative  correlations  in  Jordan  between  control  and 
exposure  sites.  Donors  with  strong  immunoreactivity  to  the 
-12  and  38  kDa  salivary  proteins  also  displayed  significant 
antibody  specificity  for  the  -30  kDa  immunodominant  pro¬ 
tein.  The  concurrent  presence  of  plasma  antibodies  to  all 
proteins  in  these  three  MW  ranges  was  the  defining  factor 
differentiating  k-means  clusters  A  from  B  (Table  4).  By 
extension,  immunoreactivity  to  all  three  salivary  proteins 
(but  not  any  one  protein  or  pair  of  proteins)  differentiated 
donors  residing  mostly  in  control  regions  from  those  donors 
residing  in  regions  with  endemic  P.  papatasi  populations.  Our 
data  showed  that  immense  antibody  specificity  to  the  -12  kDa 
protein  (but  not  the  -14  kDa  protein)  was  present  among 
many  donors.  Although  the  function  of  both  proteins  remains 
unknown,  the  differences  in  immunoreactivity  could  suggest 
divergent  functions  or  expression  patterns  or  differential  anti¬ 
body  specificity  caused  by  glycosylation.  Given  the  conserved 
pan-vector  expression  of  the  -38  kDa  salivary  protein,  the 
bloodfeeding  function  of  proteins  in  the  -30-kDa  range,  and 
the  high  degree  of  immunoreactivity  to  the  -12  kDa  protein, 
future  investigations  of  sand  fly  salivary  proteins  should  con¬ 
sider  proteins  in  these  three  ranges  in  combination. 

The  same  parameters  used  to  cluster  MENA  immunoblot 
scores  were  used  with  US  military  personnel  data.  Although 
only  two  salivary  proteins  displayed  immunogenicity  signifi¬ 
cantly  different  between  Leishmania-miected  donors  and 
unexposed  controls  (-14  and  38  kDa),  again,  antibody  speci¬ 
ficity  for  five  proteins  was  the  defining  characteristic  for  clus¬ 
ters  E-H  (Table  4  and  Supplemental  Figure  4).  Like  the 
MENA  regional  donor  data,  cumulative  immunoblot  scores 
for  the  five  proteins  (-12,  14,  26,  30,  and  38  kDa)  were  sig¬ 
nificantly  different  among  donor  members  of  each  cluster 
( P  <  0.05,  Student’s  t  test).  The  absence  of  specificity  for 
proteins  -44  kDa,  the  emergence  of  significantly  differential 
specificity  to  the  -14  kDa  protein,  and  the  apparent  immuno- 
dominance  of  the  -38  kDa  protein  (immunoreactivity  to  this 
protein  alone  defined  the  profile  of  cluster  F)  were  key  ele¬ 
ments  that  defined  the  differences  in  profiles  between  US 
military  donors  and  Egyptian  or  Jordanian  donors.  When 
assessing  the  exposure  characteristics  of  donors  in  each  of  clus¬ 
ters  E-H,  the  simple  majority  of  Leishmania- infected  donors 
made  up  the  individuals  grouped  into  cluster  F,  whereas  cluster 
G  was  populated  mostly  by  individuals  who  were  bitten  by 
sand  flies  but  not  infected  (Figure  3D).  Because  cluster  F  is 
characterized  by  the  singular  specificity  to  the  -38  kDa  salivary 
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protein,  the  average  signal  intensity  of  Leishmania-miected 
donors  in  that  cluster  to  that  protein  was  substantially  greater 
than  the  average  signal  intensity  of  the  bite  exposure-only 
donors  (2.2  versus  1.2,  respectively).  The  same  intracluster 
characteristics  were  true  between  donors  of  both  groups  in 
cluster  G.  Notably,  the  cluster  with  the  greatest  percentage  of 
unexposed  control  donors  was  cluster  H,  characterized  by  anti¬ 
body  specificity  for  ~12, 14,  and  30  kDa  salivary  proteins.  Spec¬ 
ificity  for  the  ~30  kDa  protein  occurred  in  >  50%  of  all  donors 
in  this  cluster  H,  but  among  unexposed  control  donors 
grouped  into  the  cluster,  the  split  was  an  even  50:50,  making 
immunoreactivity  to  that  protein  a  potentially  conditional 
criteria  for  cluster  definition  (Supplemental  Figure  2).  The 
dominance  of  the  -12  and  14  kDa  protein  specificity  among 
control  donors  may  suggest  that  it  is  these  proteins  that 
donors  produce  from  circulating  antibodies  as  a  result  of  expo¬ 
sure  to  other  hematophagous  arthropods  or  that  these  proteins 
exhibit  some  post-translational  modification  to  offer  epitopes 
against  which  unexposed  individuals  have  normal  circulat¬ 
ing  antibodies. 

Four  majorly  expressed  proteins,  -40,  42,  44,  and  46  kDa, 
all  matched  with  high  coverage  and  confidence  to  either 
PpSP42  or  PpSP44  (Table  3).  Both  proteins  share  84.1%  pri¬ 
mary  sequence  homology  and  resemble  yellow-related  pro¬ 
teins,  like  those  proteins  found  in  L.  longipalpis  that  function 
as  antiinflammatory  factors  on  dendritic  cells. 4S  Because  anti¬ 
body  specificities  to  both  the  -26  and  44  kDa  proteins  were 
differentiating  characteristics  between  multiple  k-means  clusters 
among  MENA  donors  (Table  4  and  Supplemental  Figure  3), 
the  expression  of  such  proteins  in  sand  flies  could  occur  in  a 
conditional  manner  to  elicit  differential  host  responses. 
Despite  that  our  data  indicate  the  -26  and  44  kDa  protein(s) 
to  be  less  immunogenic  than  the  other  three  major  salivary 
proteins  (Figure  2A  and  B  and  Table  2),  the  proposed  func¬ 
tions  and  status  in  differentiating  immunoreactivity  clusters 
make  the  case  for  additional  investigation  of  these  proteins 
as  combinatorial  diagnostic  markers  of  sand  fly  exposure. 

Enzyme-linked  immunosorbent  assays  (ELISAs)  revealed 
that  the  dominant  isotype  and  subclass  of  circulating  anti-SGS 
antibodies  is  IgG4.  Positive  correlations  were  observed  across 
k-means  clusters  for  IgGl,  -2,  and  -4  and  IgE  isotypes.  We 
also  observed  a  positive  correlation  between  concentrations 
of  IgG4  and  IgE  across  clusters  (P  =  0.02,  r  =  0.5)  (Figure  4). 
IgG4  and  IgE  are  generalized  as  markers  of  Th2  immunity 
bias  and  allergenic  type  1  hypersensitivity,  respectively.  The 
correlation  between  both  antibody  isotypes  agrees  with  prior 
reports,  where  IgE  production  is  dependent  on  Th2  responses 
and  often  associated  with  concurrent  IgG4  production.49,50 
Hence,  because  donors  displayed  circulating  anti-SGS  anti¬ 
body  profiles  with  progressive  immunoreactivity  for  all  five  of 
the  major  salivary  proteins  (Table  2),  a  possibility  emerged 
that  allergic  reactions  to  sand  fly  saliva  could  be  part  of  the 
global  host  immunological  response  dictated  by  additive 
responses  to  those  five  proteins.  These  data  further  emphasize 
the  notion  that  achieving  acquired  immunity  to  sand  fly  saliva 
may  be  the  product  of  concurrent  responses  to  multiple  sali¬ 
vary  proteins. 

To  broadly  examine  acquired  cellular  immunity,  we  moni¬ 
tored  proliferation  of  donor  leukocytes  cocultured  with  auto- 
geneic  monocyte-derived  dendritic  cells  pre-exposed  to  SGS. 
Using  a  random  subset  of  MENA  donor  samples,  we 
observed  a  trend  to  diminished  cellular  proliferation  among 


donors  from  P.  papatasi-endemic  regions  compared  with 
those  donors  residing  in  control  areas  (Figure  5A),  an  obser¬ 
vation  similar  to  the  observation  by  Rohousova  and  others.51 
Despite  the  fact  that  the  only  proliferative  responses  between 
AM  and  MA  were  significantly  different  by  Student’s  t  test 
(P  <  0.05),  the  trending  patterns  among  both  Jordanian  and 
Egyptian  donors  would  suggest  that  donors  residing  in  regions 
with  P.  papatasi  populations  elicit  a  mitigated  proliferative 
response  on  reexposure  to  sand  fly  saliva  compared  with 
unexposed  donors.  In  similar  fashion,  MDDCs  from  Iraq- 
based  military  personnel  donors  bitten  by  sand  flies  but  not 
infected  with  Leishmania  induced  significantly  less  cellular 
proliferation  compared  with  unexposed  donors  (Figure  5B). 
A  recent  report52  showed  that  CD4+  T  cells  from  sand  fly- 
exposed  humans  stimulated  with  SGS  expressed  high  levels  of 
IFN-y.  When  CD8+  T  cells  from  the  same  donors  were  stimu¬ 
lated  with  SGS,  high  amounts  of  IL-10  were  produced.  When 
both  cell  types  were  cocultured  and  SGS  was  stimulated,  only 
IL-10  was  predominantly  produced.  However,  when  IL-10 
was  blocked  during  stimulation,  cellular  proliferation  was 
substantially  increased.  The  imbalance  of  Thl  and  Th2  phe¬ 
notypes  at  the  lymphocyte  point  of  response  may  also  explain 
our  observations  in  addition  to  potential  reprogramming  of 
antigen-presenting  cells.  Moreover,  another  report  showed 
that  P.  papatasi  saliva  inhibited  the  function  of  murine  den¬ 
dritic  cells  in  mice  through  increasing  concentrations  of  free 
adenosine  nucleosides  and  AMP.53  With  the  presence  of 
the  -38  kDa  apyrase  in  our  SGS,  it  may  be  possible  that 
increased  hydrolytic  reactions  by  that  protein  to  boost  AMP 
levels  in  dendritic  cell  cultures  could  stymie  the  proliferative 
responses  of  subsequent  mixed  leukocyte  cultures.  No  signif¬ 
icant  differences  in  stimulation  indexes  were  observed  when 
donors  were  grouped  according  to  cluster  membership  (data 
not  shown).  Thus,  although  distinctive  immunoreactivity 
profiles  may  allude  to  the  biological  activity  induced  by  dif¬ 
ferential  exposure  characteristics,  those  same  profiles  do  not 
seem  to  dictate  overall  MDDC  licensing  of  immune  responses 
on  reexposure. 

Culture  supernatants  from  SGS-stimulated  US  military  per¬ 
sonnel  donors  were  used  to  assess  the  production  of  secreted 
cytokines  and  chemokines.  The  significant  difference  in  total 
cellular  proliferation  between  control  and  sand  fly  bite- 
exposed  donors  suggested  the  likelihood  of  differential  cell 
signaling  profiles  to  drive  or  inhibit  that  response.  Leishmania- 
infected  donors  produced  significantly  less  secreted  IL-8, 
IL-17,  and  IP- 10  compared  with  unexposed  control  donors 
(P  <  0.05,  Student’s  t  test)  (Figure  6).  IL-17  is  a  major  proin- 
flammatory  signaling  molecule  in  the  context  of  cutaneous 
leishmaniasis  infections  shown  to  promote  cutaneous  disease 
progression  and  severity  in  the  BALB/c  mouse  model 54  Bone 
marrow-derived  dendritic  cells  pre-exposed  to  salivary  gland 
extracts  from  the  South  American  Leishmania  sand  fly  vector 
L.  longipalpis  exhibited  diminished  IL-17  expression  when  sec¬ 
ondarily  challenged  to  inflammation.411  Moreover,  peripheral 
blood  mononuclear  cells  from  human  donors  who  self-healed 
from  cutaneous  disease  expressed  less  IL-17  compared  with 
individuals  presenting  with  active  disease  after  in  vitro  rechal¬ 
lenge  with  soluble  Leishmania  antigens.55  IL-8  and  IP-10 
(CXCL10)  are  potent  chemotactic  agents  responsible  for 
recruitment  of  neutrophils,  monocytes,  and  other  immune  cell 
types  and  have  been  shown  to  be  responsible  for  the  abroga¬ 
tion  of  the  visceral  forms  of  leishmaniasis.56  Pre-exposure  of 
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murine  cells  to  L.  intermedia  salivary  sonicates  resulted 
in  decreased  IP- 10  expression  during  subsequent  secondary 
exposures.57  The  combination  of  both  sand  fly  salivary  protein 
and  Leishmania  parasite  exposure  may  have  generated  immu¬ 
nological  memory  sufficient  to  significantly  decrease  the 
relative  expression  of  IL-17,  IL-8,  and  IP-10  after  challenge 
with  P.  papatasi  saliva  alone. 

Taken  together,  it  seems  that,  as  investigations  go  forward 
to  identify  vaccine  candidates  from  the  saliva  of  sand  flies,  a 
combinatorial  approach  will  likely  be  necessary.  Indeed, 
stronger  immune  responses  are  evident  in  hosts  pre-exposed 
to  homogenous  laboratory  colony  strains  of  P.  papatasi  than 
heterogeneous  wild  populations,58’59  with  protein  polymor¬ 
phisms  implicated  as  a  possible  explanation  for  the  phenom¬ 
ena;  this  finding  further  emphasizes  the  need  for  development 
of  a  clear  universal  baseline  profile  for  host-acquired  immu¬ 
nity  to  vector  salivary  components. 
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